Using a climate attribution statistic to inform judgments about changing fisheries sustainability.

Sustainability-maintaining catches within the historical range of socially and ecologically acceptable values-is key to fisheries success. Climate change may rapidly threaten sustainability, and recognizing these instances is important for effective climate adaptation. Here, we present one approach for evaluating changing sustainability under a changing climate. We use Bayesian regression models to compare fish population processes under historical climate norms and emerging anthropogenic extremes. To define anthropogenic extremes we use the Fraction of Attributable Risk (FAR), which estimates the proportion of risk for extreme ocean temperatures that can be attributed to human influence. We illustrate our approach with estimates of recruitment (production of young fish, a key determinant of sustainability) for two exploited fishes (Pacific cod Gadus macrocephalus and walleye pollock G. chalcogrammus) in a rapidly warming ecosystem, the Gulf of Alaska. We show that recruitment distributions for both species have shifted towards zero during anthropogenic climate extremes. Predictions based on the projected incidence of anthropogenic temperature extremes indicate that expected recruitment, and therefore fisheries sustainability, is markedly lower in the current climate than during recent decades. Using FAR to analyze changing population processes may help fisheries managers and stakeholders to recognize situations when historical sustainability expectations should be reevaluated.

Ecosystem response persists after a prolonged marine heatwave.

Some of the longest and most comprehensive marine ecosystem monitoring programs were established in the Gulf of Alaska following the environmental disaster of the Exxon Valdez oil spill over 30 years ago. These monitoring programs have been successful in assessing recovery from oil spill impacts, and their continuation decades later has now provided an unparalleled assessment of ecosystem responses to another newly emerging global threat, marine heatwaves. The 2014-2016 northeast Pacific marine heatwave (PMH) in the Gulf of Alaska was the longest lasting heatwave globally over the past decade, with some cooling, but also continued warm conditions through 2019. Our analysis of 187 time series from primary production to commercial fisheries and nearshore intertidal to offshore oceanic domains demonstrate abrupt changes across trophic levels, with many responses persisting up to at least 5 years after the onset of the heatwave. Furthermore, our suite of metrics showed novel community-level groupings relative to at least a decade prior to the heatwave. Given anticipated increases in marine heatwaves under current climate projections, it remains uncertain when or if the Gulf of Alaska ecosystem will return to a pre-PMH state.

Responses of ichthyoplankton assemblages to the recent marine heatwave and previous climate fluctuations in several Northeast Pacific marine ecosystems.

The effects of climate warming on ecosystem dynamics are widespread throughout the world's oceans. In the Northeast Pacific, large-scale climate patterns such as the El Niño/Southern Oscillation and Pacific Decadal Oscillation, and recently unprecedented warm ocean conditions from 2014 to 2016, referred to as a marine heatwave (MHW), resulted in large-scale ecosystem changes. Larval fishes quickly respond to environmental variability and are sensitive indicators of ecosystem change. Categorizing ichthyoplankton dynamics across marine ecosystem in the Northeast Pacific can help elucidate the magnitude of assemblage shifts, and whether responses are synchronous or alternatively governed by local responses to regional oceanographic conditions. We analyzed time-series data of ichthyoplankton abundances from four ecoregions in the Northeast Pacific ranging from subarctic to subtropical: the Gulf of Alaska (1981-2017), British Columbia (2001-2017), Oregon (1998-2017), and the southern California Current (1981-2017). We assessed the impact of the recent (2014-2016) MHW and how ichthyoplankton assemblages responded to past major climate perturbations since 1981 in these ecosystems. Our results indicate that the MHW caused widespread changes in the ichthyoplankton fauna along the coast of the Northeast Pacific Ocean, but impacts differed between marine ecosystems. For example, abundances for most dominant taxa were at all-time lows since the beginning of sampling in the Gulf of Alaska and British Columbia, while in Oregon and the southern California Current species richness increased as did abundances of species associated with warmer waters. Lastly, species associated with cold waters also increased in abundances close to shore in southern California during the MHW, a pattern that was distinctly different from previous El Niño events. We also found several large-scale, synchronized ichthyoplankton assemblage composition shifts during past major climate events. Current climate projections suggest that MHWs will become more intense and thus our findings can help project future changes in larval dynamics, allowing for improved ecosystem management decisions.

Effects of climate and demography on reproductive phenology of a harvested marine fish population.

Shifts in phenology are a well-documented ecological response to changes in climate, which may or may not be adaptive for a species depending on the climate sensitivity of other ecosystem processes. Furthermore, phenology may be affected by factors in addition to climate, which may accentuate or dampen climate-driven phenological responses. In this study, we investigate how climate and population demographic structure jointly affect spawning phenology of a fish species of major commercial importance: walleye pollock (Gadus chalcogrammus). We use 32 years of data from ichthyoplankton surveys to reconstruct timing of pollock reproduction in the Gulf of Alaska and find that the mean date of spawning has varied by over 3 weeks throughout the last >3 decades. Climate clearly drives variation in spawn timing, with warmer temperatures leading to an earlier and more protracted spawning period, consistent with expectations of advanced spring phenology under warming. However, the effects of temperature were nonlinear, such that additional warming above a threshold value had no additional effect on phenology. Population demographics were equally as important as temperature: An older and more age-diverse spawning stock tended to spawn earlier and over a longer duration than a younger stock. Our models suggest that demographic shifts associated with sustainable harvest rates could shift the mean spawning date 7 days later and shorten the spawning season by 9 days relative to an unfished population, independent of thermal conditions. Projections under climate change suggest that spawn timing will become more stable for walleye pollock in the future, but it is unknown what the consequences of this stabilization will be for the synchrony of first-feeding larvae with production of zooplankton prey in spring. With ongoing warming in the world's oceans, knowledge of the mechanisms underlying reproductive phenology can improve our ability to monitor and manage species under changing climate conditions.

Inferring genetic connectivity in real populations, exemplified by coastal and oceanic Atlantic cod.

Genetic data are commonly used to estimate connectivity between putative populations, but translating them to demographic dispersal rates is complicated. Theoretical equations that infer a migration rate based on the genetic estimator FST , such as Wright's equation, FST ≈ 1/(4Nem + 1), make assumptions that do not apply to most real populations. How complexities inherent to real populations affect migration was exemplified by Atlantic cod in the North Sea and Skagerrak and was examined within an age-structured model that incorporated genetic markers. Migration was determined under various scenarios by varying the number of simulated migrants until the mean simulated level of genetic differentiation matched a fixed level of genetic differentiation equal to empirical estimates. Parameters that decreased the Ne /Nt ratio (where Ne is the effective and Nt is the total population size), such as high fishing mortality and high fishing gear selectivity, increased the number of migrants required to achieve empirical levels of genetic differentiation. Higher maturity-at-age and lower selectivity increased Ne /Nt and decreased migration when genetic differentiation was fixed. Changes in natural mortality, fishing gear selectivity, and maturity-at-age within expected limits had a moderate effect on migration when genetic differentiation was held constant. Changes in population size had the greatest effect on the number of migrants to achieve fixed levels of FST , particularly when genetic differentiation was low, FST ≈ 10-3 Highly variable migration patterns, compared with constant migration, resulted in higher variance in genetic differentiation and higher extreme values. Results are compared with and provide insight into the use of theoretical equations to estimate migration among real populations.

Fine-scale population dynamics in a marine fish species inferred from dynamic state-space models.

Identifying the spatial scale of population structuring is critical for the conservation of natural populations and for drawing accurate ecological inferences. However, population studies often use spatially aggregated data to draw inferences about population trends and drivers, potentially masking ecologically relevant population sub-structure and dynamics. The goals of this study were to investigate how population dynamics models with and without spatial structure affect inferences on population trends and the identification of intrinsic drivers of population dynamics (e.g. density dependence). Specifically, we developed dynamic, age-structured, state-space models to test different hypotheses regarding the spatial structure of a population complex of coastal Atlantic cod (Gadus morhua). Data were from a 93-year survey of juvenile (age 0 and 1) cod sampled along >200 km of the Norwegian Skagerrak coast. We compared two models: one which assumes all sampled cod belong to one larger population, and a second which assumes that each fjord contains a unique population with locally determined dynamics. Using the best supported model, we then reconstructed the historical spatial and temporal dynamics of Skagerrak coastal cod. Cross-validation showed that the spatially structured model with local dynamics had better predictive ability. Furthermore, posterior predictive checks showed that a model which assumes one homogeneous population failed to capture the spatial correlation pattern present in the survey data. The spatially structured model indicated that population trends differed markedly among fjords, as did estimates of population parameters including density-dependent survival. Recent biomass was estimated to be at a near-record low all along the coast, but the finer scale model indicated that the decline occurred at different times in different regions. Warm temperatures were associated with poor recruitment, but local changes in habitat and fishing pressure may have played a role in driving local dynamics. More generally, we demonstrated how state-space models can be used to test evidence for population spatial structure based on survey time-series data. Our study shows the importance of considering spatially structured dynamics, as the inferences from such an approach can lead to a different ecological understanding of the drivers of population declines, and fundamentally different management actions to restore populations.

Stochasticity and determinism: how density-independent and density-dependent processes affect population variability.

A persistent debate in population ecology concerns the relative importance of environmental stochasticity and density dependence in determining variability in adult year-class strength, which contributes to future reproduction as well as potential yield in exploited populations. Apart from the strength of the processes, the timing of density regulation may affect how stochastic variation, for instance through climate, translates into changes in adult abundance. In this study, we develop a life-cycle model for the population dynamics of a large marine fish population, Northeast Arctic cod, to disentangle the effects of density-independent and density-dependent processes on early life-stages, and to quantify the strength of compensatory density dependence in the population. The model incorporates information from scientific surveys and commercial harvest, and dynamically links multiple effects of intrinsic and extrinsic factors on all life-stages, from eggs to spawners. Using a state-space approach we account for observation error and stochasticity in the population dynamics. Our findings highlight the importance of density-dependent survival in juveniles, indicating that this period of the life cycle largely determines the compensatory capacity of the population. Density regulation at the juvenile life-stage dampens the impact of stochastic processes operating earlier in life such as environmental impacts on the production of eggs and climate-dependent survival of larvae. The timing of stochastic versus regulatory processes thus plays a crucial role in determining variability in adult abundance. Quantifying the contribution of environmental stochasticity and compensatory mechanisms in determining population abundance is essential for assessing population responses to climate change and exploitation by humans.

Centennial-scale fluctuations and regional complexity characterize Pacific salmon population dynamics over the past five centuries.

Observational data from the past century have highlighted the importance of interdecadal modes of variability in fish population dynamics, but how these patterns of variation fit into a broader temporal and spatial context remains largely unknown. We analyzed time series of stable nitrogen isotopes from the sediments of 20 sockeye salmon nursery lakes across western Alaska to characterize temporal and spatial patterns in salmon abundance over the past ∼500 y. Although some stocks varied on interdecadal time scales (30- to 80-y cycles), centennial-scale variation, undetectable in modern-day catch records and survey data, has dominated salmon population dynamics over the past 500 y. Before 1900, variation in abundance was clearly not synchronous among stocks, and the only temporal signal common to lake sediment records from this region was the onset of commercial fishing in the late 1800s. Thus, historical changes in climate did not synchronize stock dynamics over centennial time scales, emphasizing that ecosystem complexity can produce a diversity of ecological responses to regional climate forcing. Our results show that marine fish populations may alternate between naturally driven periods of high and low abundance over time scales of decades to centuries and suggest that management models that assume time-invariant productivity or carrying capacity parameters may be poor representations of the biological reality in these systems.

A coherent signature of anthropogenic nitrogen deposition to remote watersheds of the Northern Hemisphere.

Humans have more than doubled the amount of reactive nitrogen (Nr) added to the biosphere, yet most of what is known about its accumulation and ecological effects is derived from studies of heavily populated regions. Nitrogen (N) stable isotope ratios ((15)N:(14)N) in dated sediments from 25 remote Northern Hemisphere lakes show a coherent signal of an isotopically distinct source of N to ecosystems beginning in 1895 ± 10 years (±1 standard deviation). Initial shifts in N isotope composition recorded in lake sediments coincide with anthropogenic CO(2) emissions but accelerate with widespread industrial Nr production during the past half century. Although current atmospheric Nr deposition rates in remote regions are relatively low, anthropogenic N has probably influenced watershed N budgets across the Northern Hemisphere for over a century.

Climate and population density drive changes in cod body size throughout a century on the Norwegian coast.

Understanding how populations respond to changes in climate requires long-term, high-quality datasets, which are rare for marine systems. We estimated the effects of climate warming on cod lengths and length variability using a unique 91-y time series of more than 100,000 individual juvenile cod lengths from surveys that began in 1919 along the Norwegian Skagerrak coast. Using linear mixed-effects models, we accounted for spatial population structure and the nested structure of the survey data to reveal opposite effects of spring and summer warming on juvenile cod lengths. Warm summer temperatures in the coastal Skagerrak have limited juvenile growth. In contrast, warmer springs have resulted in larger juvenile cod, with less variation in lengths within a cohort, possibly because of a temperature-driven contraction in the spring spawning period. A density-dependent reduction in length was evident only at the highest population densities in the time series, which have rarely been observed in the last decade. If temperatures rise because of global warming, nonlinearities in the opposing temperature effects suggest that negative effects of warmer summers will increasingly outweigh positive effects of warmer springs, and the coastal Skagerrak will become ill-suited for Atlantic cod.

Population diversity and the portfolio effect in an exploited species.

One of the most pervasive themes in ecology is that biological diversity stabilizes ecosystem processes and the services they provide to society, a concept that has become a common argument for biodiversity conservation. Species-rich communities are thought to produce more temporally stable ecosystem services because of the complementary or independent dynamics among species that perform similar ecosystem functions. Such variance dampening within communities is referred to as a portfolio effect and is analogous to the effects of asset diversity on the stability of financial portfolios. In ecology, these arguments have focused on the effects of species diversity on ecosystem stability but have not considered the importance of biologically relevant diversity within individual species. Current rates of population extirpation are probably at least three orders of magnitude higher than species extinction rates, so there is a pressing need to clarify how population and life history diversity affect the performance of individual species in providing important ecosystem services. Here we use five decades of data from Oncorhynchus nerka (sockeye salmon) in Bristol Bay, Alaska, to provide the first quantification of portfolio effects that derive from population and life history diversity in an important and heavily exploited species. Variability in annual Bristol Bay salmon returns is 2.2 times lower than it would be if the system consisted of a single homogenous population rather than the several hundred discrete populations it currently consists of. Furthermore, if it were a single homogeneous population, such increased variability would lead to ten times more frequent fisheries closures. Portfolio effects are also evident in watershed food webs, where they stabilize and extend predator access to salmon resources. Our results demonstrate the critical importance of maintaining population diversity for stabilizing ecosystem services and securing the economies and livelihoods that depend on them. The reliability of ecosystem services will erode faster than indicated by species loss alone.